Can competition between crystal field and Kondo effect cause non-Fermi liquid like 

behavior? 
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The recently reported unusual behavior of the static and dynamical magnetic susceptibility as well 
as the specific heat in Cei_a:Laa;Ni9Ge4 has raised the question of a possible non- Fermi liquid ground 
state in this material. We argue that for a consistent physical picture the crystal field splitting of 
two low lying magnetic doublets of the Ce 4/-shell must be taken into account. Furthermore, we 
show that for a splitting of the order of the low temperature scale T* of the system a crossover 
behavior between an SU(4) and an SU(2) Kondo effect is found. The screening of the two doublets 
occurs on different temperature scales leading to a different behavior of the magnetic susceptibility 
and the specific heat at low temperatures. The experimentally accessible temperature regime down 
to 50mK still lies in the extended crossover regime into a strong coupling Fermi-liquid fixed point. 
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PACS numbers: 71.27.+a, 71.10.Hf, 75.20.Hr, 72.15.Qm 

Introduction. The investigation of thermodynamic 
and transport properties of strongly correlated electron 
systems is of fundamental importance for our under- 
standing of elementary excitations in solid state physics. 
Especially measurements on the metallic heavy fermion 
(HF) compounds challenge the paradigm of Landau's 
Fermi-liquid concept which incorporates all lattice and 
Coulomb correlations into a renormalized quasi-particle 
mass TO* and a few Fermi-liquid parameters. 

In many cases the presence of localized moments in 
the HF compounds leads to magnetic or superconducting 
phase transitions, which either compete with each other 
in Ce based compounds or possibly even coexist as in 
Uranium based materials. The experimental evidence 
0, 01 compiled over the past ten years also indicates that 
even for HF systems with paramagnetic ground state, 
the temperature dependence of the specific heat and the 
magnetic susceptibilities often does not agree with the 
predictions of the Fermi- liquid theory, in particular when 
subject to pressure or ion substitution 4J. Therefore, 
the phenomenoiogical term non-Fermi liquid (nFL) was 
attributed to such regimes appearing in a large variety of 
different materials 

The understanding of the observed nFL behavior is one 
of the most challenging and unsolved theoretical puzzles. 
In many materials, it is ascribed to a quantum critical 
point (QCP) at which a transition temperature is sup- 
pressed to T = by an external control parameter such 
as pressure or doping^, Q. It is beheved that in the 
vicinity of such a QCP, quantum fluctuations dominate 
over thermal ones even at finite temperatures as shown 
by Hertz dl and MinisQ in a renormalized quasi-particle 
picture. Despite a tremendous experimental and theoret- 
ical effort it is, however, still not clear whether the nFL 
effects observed in HF compounds are related to novel 
low- lying non-local excitations in concentrated systems, 
true local nFL physics or simply due to competing local 



energy scales. 

Recently, the experiments showing unusual specific 
heat, magnetic susceptibility, and resistivity data for 
Cei_3;La2:NigGe4 for various concentrations have drawn 
a lot of attention since this material has the "largest 
ever recorded value of the electronic specific heat at low 
temperature" of 7(r) = AC/T « SJK^Vor^ While 
the 7 coefficient continuous to rise at the lowest experi- 
mentally accessible temperature, the magnetic suscepti- 
bility apparently tends to saturate at low temperatures. 
Experimentally, the quantum critical^ and Kondo dis- 
order scenario[3 were ruled out 01 ■ 

In this Letter, we propose a local scenario for the ob- 
served nFL behavior in Cei_a;Laa;Ni9Ge4. This is backed 
by the experimental findings that the electronic contri- 
bution to the specific heat as well as the magnetic sus- 
ceptibility normalized to the Ce concentration remains 
almost independent of the La concentration 8] . We will 
show that the competition of Kondo and crystal field 
effects leads to a crossover regime connecting incoher- 
ent spin scattering at high temperatures and a conven- 
tional strong-coupling Fermi-liquid regime at temper- 
atures much lower than the experimentally accessible 
30mK. 

Crystal-field scheme for CeA^igGe4. The Hund's rule 
ground state of Ce^+ with j — 5/2 is spHt in a tetrago- 
nal symmetry^ in three Kramers doublets. If the crystal 
electric field (CEF) parameters are close to those of cubic 
symmetry, the two low lying doublets F,.^'' and F^^-*, orig- 
inating from the splitting of the low lying Fg quartet, are 
well separated from the higher lying Fg doublet. Ignoring 
this Fg doublet, we can discuss two extreme limits. In 
a cubic environment, the CEF splitting vanishes and the 
low temperature physics is determined by an SU(4) An- 
derson model described by a strong coupling fixed point 
plus a marginal operator responsible for the particle- 
hole asymmetry jl.Qj . In a strongly tetragonally distorted 
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crystal, on the other hand, the crystal field splitting of 
the quartets is expected to be large. In this case, the low 
temperature properties are determined by an SU(2) An- 
derson model which has a significantly lower Kondo 
scale since the degeneracy N enters the denominator 
of the exponential Tk oc exp[— 1/A^J]. Then the sec- 
ond doublet at higher energies is screened at tempera- 
ture T « A = £'p(2) — i?p(i) and contributes little to the 
magnetic susceptibility. Thus, the experimental response 
would be that of a simple SU(2) Anderson model which 
was ruled out by the experiments |8l|. Therefore, we pro- 
pose that the material parameters lie in the crossover 
regime where the effective low temperature scale T* is of 
the order of the crystal field splitting A. Then, the ex- 
cited doublet will have significant weight in the ground 
state so that the total magnetic response differs from a 
simple SU(N) Anderson model. 

Formulation. Our calculation is based on an SU(N) 
Anderson model with infinite-?/ 0| whose Hamiltonian 
is given by 

H = '^Skaacl^a^kaa + '^Eaa\cta){aa-\ (1) 

ka aa 

+ 5] (l«f^)(0|cfca. + cL,|0)(afT|) , 

kaa- 

where \aa) represents a state with energy Eaa on the 
Ce 4/-shell of the a-th irreducible representation (irrep) 
with spin a and Ckaa- annihilates a conduction electron 
state with energy Skaa- transforming accordin g to the ir- 
rep a of the tetragonal magnetic point group |l^. This 
allows locally only fluctuations between an empty and a 
singly occupied Ce 4/-shell. 

We accurately solve the Hamiltonian JTJ using Wil- 
son's numerical renormalization group jlSjlUl best suited 
to deal with the competition between Kondo effect and 
CEF field splittings. The key ingredient in the NRG 
is a logarithmic discretization of the continuous bath, 
controlled by the parameter A > 1 13] . The Hamilto- 
nian is mapped onto a semi-infinite chain, where the iVth 
link represents an exponentially decreasing energy scale 
Dn A^^/^. Using this hierarchy of scales the sequence 
of finite-size Hamiltonians TCn for the A^-site chain is 
solved iteratively, truncating the high-energy states at 
each step to maintain a manageable number of states. 
The reduced basis set oi Hn thus obtained is expected 
to faithfully describe the spectrum of the full Hamilto- 
nian on the scale of Djv, corresponding to a temperature 
Tjv ^ D]^ [l^ from which all thermodynamic expecta- 
tion values are calculated. 

In order to obtain the impurity contribution to the 
specific heat, we calculate the difference between the en- 
tropy of the full model and the corresponding free 
electron gas Sfree{T), i.e. A5(r) = StotiT) - Sfree{T) 
[13. The Sommerfeld coefficient 7(r) = AC{T)/T of the 
Ce contribution to the specific heat is directly obtained 



by differentiating AS{T) with respect to T: 



AC{T) = T 



dS{T) 
dT 



(2) 



Since AC{T)/kB is a dimensionless quantity, the NRG 
predicts the absolute magnitude of AC{T) per Ce with- 
out any further fitting parameter. The experiment, how- 
ever, is needed to provide the absolute scale of the tem- 
perature axis. Since 7 has the dimensions of an inverse 
energy, we obtain an estimate of the Anderson width 
Tq — V^irp^O) by comparison with the experiment |l5l| 
{p{Q) is the density of states of the conduction electrons 
at the chemical potential). While the calculation allows 
for different values of the hybridization matrix element 
Va to incorporate different coupling strength of the mul- 
tiplets to the conduction band, we used the same matrix 
element V for all calculations to keep the number of free 
parameter as small as possible. 

Assuming a Zeeman splitting of the multiplet energy 
- o-gaP-sH, the Ce contribution to the 
magnet susceptibility is given by 



Ax 



2 2 d{S^J 
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(3) 



where the magnetic field Ha is the Zeeman splitting en- 
ergy Ha — gafJ-sH. While the g-factor is determined by 
the CEF states of the multiplets, we view them as ad- 
justable parameters and calculate Xa by applying an very 
small external magnetic field of the order oi H/Tq — 10"^ 
and estimate the values of ga by comparing with the ex- 
periment |l5| . As a consequence, the dynamical suscep- 
tibility x"('^) presented later does not contain any ad- 
justable parameter and the q integrated dynamical struc- 
ture factor S{lu) can be obtained by the dissipation fluc- 
tuation theorem, 



Siuj) [1 



MB^5aXa(w) , (4) 



where x" (w) is given by the imaginary part of the spin 
susceptibility Xa(^) =^ ^al'^a ^ i^)- Since the two 
doublets contribute differently to S{uj) for temperatures 
comparable or lower than the CEF splitting A, the pro- 
posed fit by a simple Lorentzian for x"{'^) must obviously 
fail 01 . At temperatures well above the splitting, the 
magnetic response of the two doublets is this of a quar- 
tet, and in that regime a single Lorentzian fit is justified. 
At r < A, S{u!) cannot be approximated by a Lorentzian 
consistent with the reported neutron scattering data |l(il | . 

Results. The impurity contribution to the entropy is 
plotted as function of temperature for various CEF split- 
tings A in Fig. ^ A symmetric band of D — SOFq with 
a constant density of states po = 1/2D is used in all cal- 
culations. While the entropy curves are very similar for 
small splitting, for larger A/Fq > 0.05 the decrease of 
the entropy occurs in two steps from ln(4) ln(2) and 
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FIG. 1: The impurity contribution to the entropy in units 
ks ln(2) for a ground state doublet E^(i) — — SFo and a band 

width of SOFo as function of the CEF splitting A. The inset 
shows the splitting dependent low temperature scale r*(A). 
The calculations were performed with NRG discretization pa- 
rameter A = 4, keeping A'^s — 1500 states at each NRG step. 



from ln(2) when decreasing the temperature. This 
impUes that the screening of the two doublets occurs at 
two separate energy scales once the splitting A exceeds 
the low temperature scale r*(A). We have defined the 
scale T*(A) as the temperature at which the impu- 
rity contribution to the effective spin moment A(S'^) = 
{SDtot — {SD free IS rcduccd by a factor of 2 compared to 
its local moment value of 1/4. The inset shows the split- 
ting dependency of this low temperature scale. Once 
the splitting reaches a value comparable to r*(0), the 
low temperature scale is rapidly reduced. For splittings 
slightly large than r*(0)/Fo = 0.0131, the crossover re- 
gion of the entropy from ln(4) — > is significantly ex- 
tended compared to the SU(4) curve as shown by the 
solid red (online) curve for A/Fq = 0.04. The onset oc- 
curs at higher temperature due to earlier screening of 
the upper doublet while the strong coupling Fermi-liquid 
fixed point is reached at much lower temperatures com- 
pared to a simple Kondo model. In Cei_a;Laa;Ni9Ge4 this 
temperature regime was phenomenologically attributed 
to a nFL behavior |^. Based on our calculations, we, 
however, argue that a rather extended crossover regime 
to a Fermi-liquid fixed point is observed in the experi- 
ments. 

We achieved the best agreement between theory and 
experiment in the Kondo regime of using a ground 
state doublet energy of E„{i) /Fq — —8.5 and a splitting 

of A/Fo = 0.015. We have used the comparison of the 
dimensionless experimental and theoretical specific heat 
to obtain the absolute scaling factor for the temperature 
axis, and the 7(r)-coefficient to assign an explicit value 
to Fq. Both procedures gave Fq « TIAK = 61.6meV and, 
therefore, a bare CEF splitting of A = 10.7 K was used in 
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FIG. 2: Comparison between -y{T) = C(T)/T vs T and 
the susceptibility contributions of the two doublet vs T for 
En) /To = -8.5, A/Fo = 0.015. The contribution of the 
lower doublet, Xi is much larger than the one of the upper 
doublet, X2. NRG parameters as in Fig.0 



the calculations. The corresponding entropy curve would 
be located between A/Fq = 0.01 and A/Fq = 0.04 in 
Fig.Q] No lattice renormalization effects have been taken 
into account since the experiments indicate very good 
scaling with the Ce concentration [^. The additional 
Schottky peak observed at T « 60K in the experimen- 
tal data stems from the third doublet neglected in our 
calculation. 

The comparison between the temperature dependence 
of 7(T) and x(T) is shown in Fig. |21 assuming a ratio of 
.92/51 = 2 for a good fit to the experimental data [lij |. 
The ground state doublet dominates the magnetic re- 
sponse at low temperature and tends to saturate at tem- 
peratures higher than the 7-coefficient, consistent with 
the experiments Q . We find this behavior only for CEF- 
splittings A « r*(A) while for much larger or much 
smaller values x{T) and jiT) saturate simultaneously. 

Having obtained a reasonable estimate for the ratio 
between the g values of the doublets, we can predict the 
temperature dependence of the imaginary part of the lo- 
cal magnetic susceptibility x(w). Our results are shown 
in Fig. for (a) T = and (b) at T 160 K for 
the parameter set of Fig. [21 While for high temper- 
atures (panel b) the contributions to x"{^) identi- 
cal for both doublets and can be fitted by a Lorentzian 
X"{uj) ~ Ao ^"+ff(r) of width F(T) = 13mey, this is 
not possible for O.IK < T < 40K as seen in (a). Only at 
temperatures below O.IK, Xii^) might be described by a 
Lorentzian using a F of ImeV. An effective line width of 
F — 1.2meV was used to obtain the dotted (online blue) 
line in (a) fitting x"{^)- In contrast to a naive ionic pic- 
ture, the ground state contains significant contributions 
from the first excited doublet. 

Our findings naturally explain the failure of the at- 
tempt to fit recently reported neutron scattering data 
on power samples of concentrated CeNi9Ge4 with a sim- 
ple Lorentzian for temperatures below 30K fl^. The 
change of slope in x"(w) yields a small peak for oj < 
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FIG. 3; Comparison of the contributions to x"i^) (a) T = 
AK and (b) T = 160K for the model parameters of Fig. |21 
The dotted hnes present a Lorentzian fit. Calculations were 
done with NRG discretization A — 2.5 keeping Ns = 6400 
states at each step. 
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FIG. 4: The g-integrated dynamical structure factor S{ui) for 
different temperatures versus frequency for NRG parameters 
as in Fig. The inset shows an estimate for temperature 
dependent relaxation rate r(r) fitted to the lower doublet 
only. 



in S{u;,T = AK) shown in Fig. H Note, that the un- 
usual behavior of Xa{'-^) in the crossover regime can 
also not be explained as originating from a distribu- 
tion of Lorentzians [13 used to fit S{uj) in [3. Such 
phenomenological approaches contain obviously limited 
amount of information if no consistent physical picture 
for all physical properties emerges. A rough estimate 
of r(r) plotted as inset to Fig. ^ indicates F cx T as 
in the experiment 0| before it saturates below lOOmK. 
Our absolute values for F(T) are roughly a factor of two 
larger than those reported in . This is consistent with 
the fact that T* for Ceo.5Lao.5Ni9Ge4, the material we 
used to fix our model parameters, is larger that T* for 
CeNigGe4. In addition, our error in the estimate of the 
CEF splitting A sensitivly determines the low temper- 



atur scale T* as depicted in the inset of Fig. ^and there- 
fore our absolute energy scale. 

Summary and discussion. A simple physical picture 
of the competition between Kondo effect and CEF split- 
ting leads to an extended crossover region from the high 
temperature free multiplet to the low temperature strong 
coupling fixed point if the crystal field splitting is of the 
order of the low temperature scale T* . Both doublets 
contribute significantly to the magnetic response at the 
fixed point yielding different contributions to the static 
and dynamic susceptibility. We propose that the origin 
of the nFL behavior in Cei_a;Laa;Ni9Ge4 is related to this 
extended crossover region compared to a simple Kondo 
model or a two channel Kondo lattice scenario 0, . 
This provides a consistant picture for the temperature 
dependence of specific heat and the magnetic response in 
agreement with all experimental data in all regimes. We 
hope to include lattice coherence effect in future calcula- 
tions to explain the transport properties as well. 
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